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The field theory of a fluctuating d-wave superconductor is constructed and proposed as an effective 
description of superconducting cuprates at low energies. The theory is used to resolve a puzzle posed 
by recent experiments on superfluid density in severely underdoped Y Ba2Cu30&+ x - In particular, 
the overall temperature dependence of the superfluid density at low dopings is argued to be described 
well by the strongly anisotropic weakly interacting Bose gas, and thus approximately linear in 
temperature with an almost doping-independent slope. 



The superconducting state of underdoped high- 
temperature superconductors is long known to be anoma- 
lous. Whereas the pseudogap temperature T* is high, 
and only increasing with underdoping, the superconduct- 
ing transition temperature T c , together with the zero- 
temperature superfluid density p(0) at the same time is 
continuously vanishing This is in stark contrast to 
the predictions of the standard BCS theory, in which of 
course T c and T* would be essentially identical, and p(0) 
proportional to the electron density, i. e. p(0) *~ 1 — a;, 
where x is the number of holes per lattice site (doping). 
It has recently been argued , |3 , that when pairing in 
the d-wave channel is accompanied by strong repulsion 
the amplitude of the order parameter |A| and p(0) may 
behave oppositely as half-filling (x = 0) is approached, 
with the divorce of T* and T c arising as a natural conse- 
quence. Disordering of the d-wave superconductor (dSC) 
by quantum fluctuations thus appears, at least in this 
respect, to be similar to doping of holes into the Mott in- 
sulator (MI) as described by the effective gauge theories 
of the t-J model Q. 

There exist, however, experimental results that ap- 
pear to contradict the theories of the fluctuating dSC 
13 > @- 0- or the gauge theories of the t-J model Q. 
Whereas these approaches correctly yield p(0) ~ x as 
has been seen, they also generally imply that the 'ef- 
fective charge' of quasiparticles becomes proportional to 
doping and thus dp(T) / cIT\t->o ~ x 2 , which has not. 
On the contrary, the experiments in strongly underdoped 
single crystals and thin films ^(j of YBo,2Cu30q+ x 
(YBCO) show the superfluid density to be approximately 
linear in temperature over most of the temperature range, 
with an almost doping-independent slope, and with a 
higher power-law emerging near T = 0. The critical re- 
gion, almost 10 Kelvins wide at optimal doping, appears 
also to have completely disappeared at low dopings. The 
doping-independent slope becomes particularly troubling 
in light of the recent measurements of heat transport in 
high-purity single crystals of YBCO 11], in which gap- 
less quasiparticles in strongly underdoped samples indeed 
appear decoupled from the external magnetic field. This 
would agree with the theoretical prediction that quasi- 
particles should gradually become 'neutral' with under- 
doping, but at the same time makes the behavior of the 
superfluid density all the more puzzling. 

In this Letter I present the field theory of the quantum- 



fluctuating d-wave superconductor (dSC). The spin sec- 
tor, which was the primary focus of the earlier studies |(|, 
0j E3> nas been embedded into a more general frame- 
work describing both spin and charge degrees of freedom 
at finite dopings. The theory describes nodal quasipar- 
ticles coupled to the phase fluctuations of the order pa- 
rameter. Quantum fluctuations should arise from the 
Coulomb interaction, which becomes increasingly detri- 
mental to phase coherence near half-filling 3]. Such an 
effective theory should provide a correct description at 
energies below the amplitude |A|. From this starting 
point two main results are further derived. First, the 
general theory is shown to be considerably simplified by 
being rewritten in terms of the dual variables. Remark- 
ably, in this form it resembles the effective SU(2) gauge 
theory of the t-J model 0. Second, the transformed 
theory facilitates a simple, and a physically transpar- 
ent calculation of the ab-plane superfluid density p(T) 
in the relevant doping regime. Although the quasiparti- 
cle 'charge' is found to indeed be gradually vanishing, on 
the scale of T c the region of temperatures where the slope 
of p(T) would become too strongly doping-dependent be- 
comes negligible at low dopings. In the leading approxi- 
mation, the overall form of p(T) is given by the conden- 
sate of the strongly anisotropic weakly-interacting three- 
dimensional Bose gas. This, on the other hand, is shown 
to be Ap(T) ~ Tln(T/t), with t as the Josephson inter- 
layer coupling, and thus also approximately linear in tem- 
perature, but with the slope at most logarithmically de- 
pendent on doping. This also explains the absence of a 
discernible critical region and a slight curvature of the 
data near T — (Fig. 1). Finally, properties of the 
insulating phase at x = are briefly discussed. 

Let me begin by postulating the action for the low- 
energy quasiparticles of a two-dimensional (2D) phase- 
fluctuating dSC, S = J^ T dr J d 2 xL, with the La- 
grangian density L — + L$, and 

Lis = *i[7o(<9 r - ia o) + VFli{d x - ia x ) + (1) 
VA72{d y - ia y )]^i + (1 — > 2, x <-> y) + iJ^v^ + 

£• = ^ p (a^v^(A-At) T + ^( V ^ + A,) 2 (2) 
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Two four-component fermionic fields ^1,2 describe the 
gapless spin- 1/2 excitations near the two pairs of di- 
agonally opposed nodes, vp and v& ~ |A| are the 
two characteristic velocities of the low-energy spectrum, 
{7^7 lv] = 28p V , and and ^nl^n, H = 0, 1, 2 are the 
charge and spin currents, respectively, as defined in 0- 
Whereas the Fermi velocity vp should be approximately 
independent of x, va may be assumed to be a decreasing 
function of doping Q- a p and are the Berry and 
the Doppler U(l) fields that furnish the coupling of the 
quasiparticles to the fluctuating phase of the order pa- 
rameter ,6] ; the integration over the auxiliary gauge fields 
constrains €p l/ pd v (a p ,v p ) = 7r(J $1 -J$ 2 , J $1 + J$ 2 )p, 
where J^> 1 2 are the vortex current densities. Kq and 
K12 = K ~ Ep ~ 1 — x are the bare compressibility 
and the bare superfluid density, respectively, which de- 
rive from the integration over the high-energy fermions. 
Ap is the physical electromagnetic vector potential, and 
h the bare chemical potential. The parameter a tunes 
quantum fluctuations, and /3i t 2 > describe the short- 
range repulsive interactions between vortex loops. Terms 
that are irrelevant at low-energies ^| have been omitted. 

At h — L$ represents the continuum limit of the lat- 
tice theory discussed in 0. Its form may be also under- 
stood on the basis of symmetry: 1) the usual electromag- 
netic gauge invariance under Ap — > Ap + dpX, 2) the in- 
ternal gauge invariance under a M — ► ap + dpX, 3) the Ising 
symmetry under $1 <-> $2, 0,p — o,^, spin up<->down, 
and 4) the gauge invariance under A^ — > A^ + dpX ± , 
together with the requirement of analyticity in <& 1)2 dic- 
tate the form of L as unique to the lowest non-trivial 
order in the fields and their derivatives. An important 
feature of the Lagrangian is the addition of the chemi- 
cal potential h, which is necessary to allow for a finite 
doping. It is introduced by shifting Ao —> Aq + ih, af- 
ter which it has been absorbed into a redefined Vq, as 
vq + ih — > vq. It then appears only as the fictitious exter- 
nal 'magnetic field' in the r-direction in L$. As will be 
discussed shortly, having h ^ is crucial for obtaining 
the correct charge dynamics of the fluctuating dSC. 

The form of L may also be justified on phenomenologi- 
cal grounds, since, as discussed below, it describes a novel 
Ml-dSC transition of possible relevance to cuprates. It is 
convenient, however, to derive the equivalent dual form 
of L$, more amenable at h 7^ 0, first. Duality is most 
precisely established on a lattice, and for the 'hard-spin' 
version of the complex fields. Using a fairly standard set 
of transformations Q it is easy to show that, modulo 
analytic terms, 
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where r labels the sites of the 2+lD quadratic lattice, 
and A and Ax are the lattice gradient and the curl. 
Taking the continuum limit and going into the 'soft-spin' 
representation 0] this implies that 



L* = 4rK + A^) 2 + hJ2 K(do - i(vo + (-)"ao))M4) 
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L$ may therefore be alternatively understood as describ- 
ing the coupling of the fields a p and v p to two non- 
relativistic bosonic fields of unit electromagnetic charge, 
which are dual to the original vortex fields $1,2- It is 
interesting to note that, although presumably more gen- 
eral, Eq. 4 is similar to the conjectured effective theory 
of the s-flux phase within the SU(2) gauge theory of the 
t- J model |4| . Somewhat similar connection between the 
U(l) gauge theory of the t-J model and the dSC with 
topologicallytrivial phase fluctuations 5] has been noted 
in the past [lj|, jlfij - 

Consider the superconducting phase described by L. 
Assuming a > 0, /3i.2 > and minimizing for h > 



2a one finds |(6i)| 2 = |(6 2 )| 2 = (h 2 - h 2 c )/2(2f3 1 + fa). 
To the quadratic order L$ then reduces as 



- - + A p ) 2 + 9 -^{vl + a 2 ) - ihn b v , (5) 
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where Pb(T) and nb are the superfluid and the total 
density of the bosons, respectively. Since both U(l) 
fields are massive, nodal quasiparticles interact only via 
short-range interactions, and therefore represent well- 
defined low-energy excitations, conducting heat, for ex- 
ample. At low energies L = L sp + L c h, where the spin 
part of the Lagrangian to the leading order is L sp = 
Lis - iJp(vp + Ap) + (p b /2)a 2 p [H, and 

Lch = iJp(v+A)p + ^(v+A)l + ^vl~ i hn b v . (6) 

Setting Aq = and integrating over vq one finds 

hnb 



T 2 



2(K + p b (T)) K +p b (T) 
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iJ.{v + A) + ^[v + Af + ^l^, 



where X = (Xi,X2). The first term describes a short- 
range repulsion between fermions and as such is irrelevant 
at low energies. The second determines the renormalized 
chemical potential: p{T) = —hnb/(Ka + pb(T)). Since 
x ~ — /i(0), one finds nb ~ x. The rest of L c h determines 
the superfluid density. Integrating v generates the term 
zJ ■ A in L c h, with the ratio z = pb(T)/(K + pb(T)) as 
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FIG. 1: H c i in an underdoped YBCO J3 converted 
to ab-plane superfluid density in Kelvins |17|1 as p = 
(4irfr cl /((ln(K)+Q.5)20.7Gouss)K6.21iQ(c/10A), for k = 100 
and the c-axis anisotropy 10~ 4 |20| . The line is pb(T) from 
the Eq. 9. with p b (0) = 28.3JC. 

the Fermi liquid 'charge renormalization' parameter. The 
electromagnetic field gradually decouples from quasipar- 
ticles with underdoping, in accord with the observed in- 
sensitivity of thermal conductivity to the magnetic field 
in underdoped YBCO [n). 

The suppression of the quasiparticle charge, and par- 
ticularly its effect on the superfluid density, may be 
alternatively understood by integrating fermions in L 
before v. At T ^ 0, this reduces K in Eq. 7 as: 
K -► KIT) = K - (2hx(2)/n)(v F /v A )T + 0(T 2 ). In- 
tegrating then v yields the total superfluid density: 

p{T)- l =K{T)- l + Pb (T)-\ (8) 

The last result, also known as Ioffe-Larkin rule p|, is eas- 
ily seen to hold both in the underdoped (large K), and 
in the overdoped regime (large Pb(0)). At T = 0, there- 
fore, p(0) — Kp b (0)/(K + jOft(O)) ~ n b ~ x at low x, 
which agrees qualitatively with more microscopic calcu- 
lations Q, as we H as with experiment At T 0, 
however, p b (0) - p b (T) cx T 3 in 2D [HI, and the leading 
temperature dependence of p(T) at low T comes from 
the first, quasiparticle, term in Eq. 8. Differentiation 
gives dp(T)/dT = (p(0)/K) 2 dK(T)/dT ~ z 2 at low T. 
Such a strong, ~ p 2 (0), dependence of the slope with 
doping, however, appears to be clearly contradicting the 
experiment 9], which defines our main problem. 

While the reasoning from the previous paragraph is 
correct in principle, the range of temperatures where it 
applies clearly depends on the specific form of the bosonic 
term in Eq. 8. One needs therefore to understand p b (T) 
in the underdoped region, which corresponds to the di- 
lute system of bosons, in more detail. Before turning to 
the case of cuprates, let me first briefly review a canoni- 
cal example of bosons with a mass m interacting with a 
two-body interaction V{x) = \5{x). Dimensional anal- 
ysis dictates that the superfluid density in an isotropic 
3D system may be written as p b — r~ 3 G(m\/r,mTr 2 ), 
where r is the average distance between particles, and 



G(x, y) a dimensionless function of two dimensionless ar- 
guments, with G(0,0) = 1. The dilute limit r — > oo 
is equivalent therefore to the weakly interacting high- 
temperature limit, where G pa 1 - (T/T BEC ) 3/2 . This 
is just the temperature dependence of the condensate 
of non-interacting bosons, with Tbec being the Bose- 
Einstein condensation temperature. At low tempera- 
tures, of course, interactions always turn the temperature 
behavior into ~ T 4 (in 3D), but this higher power-law 
sets in only below the characteristic interaction tempera- 
ture scale AT\ ~ {m\/r)TBEC- Similarly, the critical re- 
gion AT C - {mX/r) 2 T BEC , so both AT A , C < T c w T BEC 
in the dilute regime. This is all just another way of phras- 
ing the irrelevancy of short-range interactions near the 
quantum critical point of 3D bosons. Short-range inter- 
actions are irrelevant in 2D as well but only marginally 
so, and consequently the weakly interacting regime is 
harder to reach. Nevertheless, in either case the bosonic 
superfluid density in the sufficiently dilute regime ap- 
proaches the characteristic temperature dependence of 
the condensate, Ap b (T) ~ T D I 2 , except in the narrow 
low-temperature and critical regions |l8j . 

What is the function p b (T) then in severely under- 
doped cuprates? First, the phase coherence in the su- 
perconducting state is known to be strongly anisotropic, 
but nevertheless fully three-dimensional. We may in- 
clude this feature by additionally coupling the different 
superconducting layers, each described by Eq. 4 and la- 
beled by I, via a weak Josephson coupling of the form 

'Ei, n =y K,ib n ,i+ii witn 1 < Pb(0)- Second, in such 
a strongly anisotropic system the Coulomb interaction 
within a layer is screened well by the other layers, and 
except at extremely low densities p b {0) ~ t it may in- 
deed be taken to be effectively short-ranged [19|. as in 
Eq. 4. Away from the ultimate quantum critical region 
of the theory (1) where both the long-range nature of the 
Coulomb repulsion, and possibly the gauge field fluctua- 
tions, would finally become important, such a short-range 
interaction is then irrelevant, and Pb{T) to the zeroth 
order in interaction is just the condensate of the non- 
interacting system. In the anisotropic 3D case this gives 
p b (T) = Pb (0) - TF{t/T)/(2w), where 

F(x) = - dzhx(l-e- xsin2 (™W), (9) 
Jo 

with F ps 2.6f2(T/7rf) 1 / 2 for T«f, and F pa \n(T/t) 
for t < T, Since t < /0&(0), on the scale of T c p b (T) 
behaves approximately linearly with temperature. This 
being the case, it is now the bosonic term in Eq. 4 
that determines the form of p(T) everywhere, except at 
temperatures below ~ tx 4 . Furthermore, assuming the 
T = c-axis superfluid response to be also determined 
by the bosonic component fixes the parameter t: at low 
dopings t = 2abp c (0)/(xc 2 ), where x/(abc) is the bo- 
son density, a, b and c are the dimensions of the unit 
cell in the a-, b-, and c-directions, and p c (0) the c-axis 
superfluid density (in Kelvins 0])- As an illustration, 
taking x pa ab/c 2 pa 0.1 and the measured value [2(j 
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of p c (Q)/p(0) = 10 -4 , a satisfactory fit to experimental 
data [9| is achieved even by neglecting the quasiparticle 
term in Eq. 8 completely, and by adjusting only pb(0) 
(Fig. 1). In the overdoped regime, on the other hand, 
Pb(T) » K(T), and p(T) should cross to the standard 
BCS result, as may have been already observed plj . 

For completeness, I also briefly describe the non- 
superconducting phase of L. For h < V2a, = (62) = 
0, and the bosons are in the incompressible phase. The 
integration over the bosons then yields 



L<T> 



^i(^ + A M ) 2 + (W^p) 2 + (w^p) 2 ) (|()) 
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to quadratic order, where m 2 ~ a + Oifix^). At low 
energies one may still write L = L sp + L^, but now 
with L sp as the three dimensional quantum electrody- 
namics (QED 3 ) @, 0- Quasiparticles interact via long- 
range gauge interaction and cease to be sharp excitations. 
Furthermore, fermions may acquire a very small mass, 
< 10~ 4 mb |22j . which would imply the antiferromagnetic 
ordering in the system. L c h, on the other hand, after the 
integration over becomes 
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The first term implies that the system is an insulator 
with a charge gap ~ m^. The second term after a partial 
integration may be rewritten as ~ Jp,€pupd^B p , where 
Bp = e P p V dpA Vl so a uniform magnetic field becomes 



completely decoupled in the insulator. Since the life- 
time of fermions is inversely proportional to the fcrmion 
mass, and thus likely to be very long, thermal conductiv- 
ity in the insulator should essentially remain the same. 
This could explain why the thermal conductivity of even 
a weakly insulating YBCO apparently remains linear in 
temperature [Til ]. 

It may be noteworthy that L also has a metastable state 
for h > he, with (61) = 0, and 2|(6 2 )| 2 = (h 2 - /i 2 )/(/3i + 
/3 2 ). In this state L$ — * pb(vp + a M ) 2 /2 — ihn b (v Q + a ), 
to quadratic order. Although pb 7^ 0, the full system is 
then actually metallic, since fermions acquire back their 
full electromagnetic charge and form four hole pockets 
with a small Fermi surface, ~ h. 

To summarize, the proposed effective theory predicts 
that a weakly fluctuating d-wave BCS superconductor 
at large dopings crosses over to the strongly phase- 
fluctuating regime in the pseudogap, low-doping regime, 
and ultimately suffers a transition into the insulator with 
possible antiferromagnetic order. Superfluid density as 
a function of temperature at low dopings becomes well 
approximated by the Bose condensate in the strongly 
anisotropic 3D Bose gas, and thus appears linear with 
an essentially doping-independent slope over most of the 
temperature range. 
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